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Abstract 

Mandatory emission trading schemes are being established around the world. 
f-*) ' Participants of such market schemes are always exposed to risks. This leads 

CN ■ to the creation of an accompanying market for emission- linked derivatives. To 

^ i evaluate the fair prices of such financial products, one needs appropriate models 

^ ^ I for the evolution of the underlying assets, emission allowance certificates. In this 

paper, we discuss continuous time diffusion and jump-diffusion models, the latter 
. enabling one to model information shocks that cause jumps in allowance prices. 

We show that the resulting martingale dynamics can be described in terms of 
■ non-linear partial differential and integro-differential equations and use a finite 

p I . difference method to investigate numerical properties of their discretizations. 

^ I The results are illustrated by a small numerical study. 
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; 1 Introduction 

m 

Emission Trading Schemes (ETSs) have recently been designed with the main aim to 
O I reduce emissions of greenhouse gases and other pollutants. Two examples of ETSs are 

the EU ETS (European Union Emission Trading Scheme) and the US Sulfur Dioxide 
Trading System. In such schemes, the regulator allocates each market participant a 
number of credits, each of which gives the holder the right to emit a specified amount 
of pollutant (e.g. CO2). At the end of a compliance period, each agent should not 
have released more pollutant than covered by the credits it holds at that time, or it 
will have to pay a fine proportional to the excess of the pollutant amount over the one 
corresponding to the credits held. During the compliance period, however, each agent 
can trade credits with other market participants, depending on whether it is cheaper 
to reduce emission or to buy credits. 

The present study relies on a detailed mathematical model for such a scheme based 
on market equilibrium. We exploit its implications for the risk neutral allowance price 
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evolution. The key issue is a feedback relationship between allowance prices and pol- 
lution abatement activity. Namely, any increase in allowance price enforces emission 
reduction since agents would then tend to sell their allowances. Hence the increase of 
the allowance price encourages supply of certificates and lowers the probability of non- 
compliance, which, in its turn, tends to diminish allowance prices. Following [15], we 
describe this feedback relationship in terms of market fundamentals, using a stylized 
evolution of the expected non-compliance and emission savings costs. We pay partic- 
ular attention to discontinuities in the information flow, assuming that some events 
may cause market participants to change their views on the future reduction volumes 
required to reach compliance. We show how risk neutral allowance price reacts to such 
"information discontinuity" and study this effect in derivatives valuation. 

The paper is organized as follows. In Section |2], we extend the discrete-time frame- 
work from [15] to the continuous-time setting. In Section |3l we interpret our results 
in the context of diffusion processes, derive the corresponding partial differential equa- 
tions and discuss option pricing. Section H] is devoted to modeling shock events in the 
information flow. Our description is based on jump-diffusion processes and requires 
solving a nonlinear partial integro-differential equation. We show that solutions to 
these non-linear equations from Section 3 and 4 satisfy the maximum principle. In 
the last section we propose a numerical implementation of its solution. In particular, 
we prove that the discretized problem possesses a unique solution and satisfies the 
maximum principle. 

The suggested allowance price model is suitable for option pricing. It turns out 
that the fair price dynamics of European options written on emission allowance prices 
can be obtained by solving a linear partial integro-differential equation. We illustrate 
our results by numerical examples. 

The literature on this subject is rather extensive, and we refer the interested reader 
to a nice expository work [27] which reviews the fundamental concepts of the environ- 
mental economics and provides a valuable guide to publications, which, however, is 
far from being complete. The economic theory of allowance trading goes back to [12] 
and [21] , where the authors proposed trading the public good environment by means 
of transferable permits. Important results in dynamic allowance trading were obtained 
in [11], [28], [23], [19], [23], [26], [20] and in the literature cited therein. Recently, after 
the introduction of the real-world emission market EU ETS, the empirical evidence 
has become available. The experience gained from the operation of Phase I of the EU 
ETS is discussed in [13], and a detailed analysis of spot and futures allowance prices 
from this market is given in [29] and [30]. The contributions [2] and [22] are devoted to 
econometric modeling of emission allowance prices. Beyond confirming stylized facts of 
financial time series for prices of emission allowances, Markov switch and AR-GARCH 
models are suggested. The modeling of dynamic price equilibrium is addressed in [3] 
and [5], which provide a mathematical analysis of the market equilibrium and use opti- 
mal stochastic control theory to show social optimality of emission trading schemes. A 
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recent work considers equilibrium of risk averse market players and elaborates on 
risk neutral dynamics. The problems of derivatives valuation in emission markets are 
also addressed. The paper [8] discusses an endogenous emission permit price dynam- 
ics within equilibrium setting and elaborates on the valuation of European options on 
emission allowances. The dissertation [31] and the paper [25] deal with the risk-neutral 
allowance price formation within the EU ETS. In that paper, when utilizing equilib- 
rium properties, the price evolution is treated in terms of marginal abatement costs 
and optimal stochastic control. The work [7] is also devoted to option pricing within 
EU ETS. The authors suppose that the drift of allowance spot prices is related to a 
hidden variable, which describes the overall market position in allowance contracts, and 
make use of filtering techniques to derive option price formulas which reflect specific al- 
lowance banking regulations valid in the EU ETS. Finally, the recent work [B] presents 
an approach where emission certificate futures are modeled in terms of a deterministic 
time change applied to a certain class of interval-valued diffusion processes. 



2 Modeling emission markets in discrete time 

During the compliance period (typically 3-5 years), each ETS participant dynamically 
adjusts its production processes (and hence its emissions) and trades emission credits 
at times t = 0,1,2, ... ,T < oo to maximize its revenue. In this setting, allowance 
price reaches its equilibrium determined by supply and demand of emission allowances. 
In what follows, we base on the model from [13] which characterizes the equilibrium 
allowance prices in terms of non-compliance uncertainty and abatement costs. This 
characterization forms the starting point for our analysis. 

Let (yi, J^,P , {J^t}t=o) be a filtered probability space. We assume that all the 
processes considered in this section are adapted to {J^t}t=Q- Consider a market with 
a finite set / of agents who must comply with the ETS rules. Assume that, for each 
2 G /, an exogenously given stochastic process {EI}JSq describes the so-called "business 
as usual" emission of agent i. That is. El = El{uj) stands for the total pollution of 
the agent i which will be emitted during the time interval (t, t + 1] if no abatement 
measures are applied by the agent. Suppose that each agent i can decide, at any 
time t = 0,...,T — l,to perform a reduction of Q pollutant units to be emitted during 
{t, t+V\. The cost of abatement is modeled by a function of the reduced volume that can 
be random (reflecting the uncertainty in fuel prices). Thus, if at time t = 0, . . . , T — 1 
agent i decides on reduction by x G [0, oo) units for the time interval (t, t + 1], then it 
costs it Cl{x) = Cl{x){u), where Q : [0, oo) x 1] t-)- [0, oo) is B{[0, oo)) measurable, 
and for each w G fi, the mapping x t— )■ Cl{x){uj) is strictly convex and continuous with 
Ct(0) = 0. Therefore, following this abatement policy {Q}J=q, agent i will have born 
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by the compliance date T the total abatement costs of 

T-l 



(1) 



t=0 

For each i&I,u&Q, t = 0,...,T — 1 and a G [0, oo), we introduce the "reduction 
volume" 

rj(a) = ri{a){uj) := argmax [ax - Ci{x){uj) : x G [0, EKu)]} . (2) 

This quantity gives the "locally optimal" reduction volume for agent i for the time 
period (t, t + 1] given that the price of one allowance unit is equal to a for that time 
period; we may assume (and Proposition 12.11 (b) below supports this) that, being 
rational, agent i will implement emission reduction at that level. As is well known (see 
e.g. [IS]), under the above assumptions rl{a){uj) is a non- decreasing and continuous 
function of a G [0, oo) for each w G and t = 0, . . . , T — 1. Next denote the total t-th 
time period reduction by 

rt(a) := J]rj(a), aG[0,oo), (3) 

which represents the total reduction by all agents in the market for the time period 
(t, t + 1], given the time t price of one allowance unit is a. 

Suppose that, at any time t = 0, . . . ,T, credits can be traded at the spot price Af. 
Denote by 'd] the change in the allowance number held by agent i at time t. Then, 
given the allowance prices {At}J^Q, the position changes {i!^l}f=Q will result for agent i 
in the total trading costs of 

T 

Y.^IA. (4) 

t=o 

Here and in what follows, for simplicity's sake we assume, zero interest rates (or, 
equivalently, that all the prices are already discounted). 

Further, the total pollution of agent i during the compliance interval [0, T] can be 
expressed as the cumulative business-as-usual emission less the agent's total reduction: 

T-l T-l 
t=o t=0 

Denoting by 7i the agent's initial allowance allocation, we observe that i will hold 
7j + ^21=0 allowances at time T and hence its loss resulting from the potential 
penalty payment at rate vr (the penalty for the amount of emissions corresponding 
to one allowance) is 



"T-l 

TT 

.t=0 



J2iEl-Ct-^l)-f-^'T ■ (5) 
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Finally, we define the space of feasible trading strategies i?* = and abatement 

strategies = of agent i E I as 

U':={{^\C) :0<Q<EI t = 0,...,T-l}. 

In view of ([1]), and ([5]), the total revenue of agent i following an admissible policy 
e is equal to 

T-l VT-l 1 + 

L^.X^\ C) := - Y.(^lAt + C\Q) - ^'^At - n ^(E^ - Q - d\) - 7' - . 

t=0 Li=0 

To specify risk preferences, we describe agents' risk attitudes by individual utility 
functions U\ i E I, that are assumed to be continuous strictly increasing and concave. 
For a random variable X, consider the utility functional 

X^u'{X) = E{U'{X)), 

which is defined whenever the expectation exists and is not —00. Given an allowance 
price process A = {At}J^Q, agent i behaves rationally in the sense that it maximizes 
its utility from the terminal wealth 

by an appropriate choice of the strategy that we denote by (^?**,,^**). Following the 
standard apprehension, a realistic market state is described by the so-called equilibrium 
— a situation where the allowance price, positions and abatement measures are such 
that each agent is satisfied with the own strategy and, at the same time, the natural 
restrictions are met. 

Definition 1. An adapted process A* = {A*}f^Q is called an equilibrium allowance 
price process if, for each i E I, there is a strategy ("i?**,^**) E W such that we have 
u\L^'^\'d*\C')) < 00 and 

(i) the cumulative changes in positions are in zero net supply: 

"^^f = Qfor allt = 0,...,T, 

(ii) each agent i E I is satisfied with its own strategy in the sense that, for each 

^ such that ^*)) exists, one has 
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In [15], this equilibrium notion was used to establish a reduced-form model which 
describes the allowance price evolution from the risk-neutral perspective. This ap- 
proach utilizes the following three properties of the above equilibrium. 

(a) There is no arbitrage since any profitable strategy would immediately be followed 
by all agents. 

(b) Given a technology with lower reduction costs than the present allowance price, 
it is optimal to immediately reduce one's pollution and take profit from selling 
allowances. 

(c) There are only two final outcomes for allowance price: at maturity, either the price 
will vanish if there is an excess in allowances, or, in the case of their shortage, 
the price will rise to the penalty level. Exact coincidences of allowance demand 
and supply at maturity occur with zero probability under broad assumptions and 
can be neglected. 

The last property follows from the assumption that the random variable Yli 
representing the total emissions at maturity has a continuous distribution, given the 
information up to time T — 1. Under mild additional assumptions, [15] claims that the 
above assertions can be deduced from the equilibrium in the following form. 

Proposition 2.1. Suppose that {A^}J^q is an equilibrium allowance price process and 
{Q^}I=o J i G /, are corresponding equilibrium abatement strategies. 

(a) There exists a measure Q on {Q, J-") which is equivalent to P and such that 
{A^}J^Q is a Q-martingale. 

(b) For each i E I one has 

Ct' = r\{Al), t = 0,...,T-l, (6) 
with the reduction functions rl, t = 0, . . . ,T — 1, from ([2]). 

(c) The terminal value of the allowance price is given by 

= ^1 (y. - ^r) - Y j > j . (7) 

where l(-B) is the indicator of the event B. 

In fact. Proposition 12.11 states the above-mentioned feedback relationship. Namely, 
the equilibrium allowance price process ^ martingale under Q ~ P that 

has the terminal value ([7]). However, this terminal random variable depends on the 
intermediate values through ([6]). A surprising consequence of the feedback 
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relationship is that, from the risk-neutral perspective, only the cumulative market 
quantities are relevant. To see this, introduce the overall "business- as-usual" allowance 
shortage by 

iei \t=o 

Further, recall from (|3]) the cumulative abatement function rt{a) to express the risk- 
neutral certificate price dynamics using ([7]) and the martingale property of {AI}]Lq 
under Q as 

T-l 



a: 



TT 



t=0 



t = 0, 



.T-l. 



Although the individual market attributes seem to be irrelevant here, the reader should 
notice that this picture appears only from the risk- neutral viewpoint. 

With this, the problem of risk neutral allowance price modeling boils down to the 
following task: 



Given a measure Q ~ P and reduction functions 
{rt}J~Q, describe the random variable £t and 
determine a Q-martingale {A^}J^q such that 

A*^ = nl [St - Eto rMl) > 



In this form, transition to the continuous time case is straightforward and the resulting 
problem can be stated as follows: 



Given a measure Q ~ P and reduction functions 
{''^t}te[Q,T], describe the random variable £t and 
determine a Q-martingale {A^}t^[Q^T] such that 

A*T = TTl (St - J^rs{A:)ds > o). 



Problem (jH]) is the starting point of our investigation in this paper. This approach 
utilizes the ingredients {rt}te[o,T] and St, which is reasonable from the practical per- 
spective since the price-dependent reduction functions {rt}te[o,T] can be estimated from 
the market data and the potential allowance shortage St can be modeled in terms of 
emission fluctuations. 



3 Continuous time case: diffusion models 

Modeling in continuous time stipulates that the compliance period is an interval [0,T] 
and that all the relevant random evolutions are described by adapted stochastic pro- 
cesses on a filtered probability space {Q, J^,Q, {J^t}t£[o,T]) which is equipped with a 
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"spot martingale" probability measure Q ~ P. Given a random variable 8t and ap- 
propriate non-decreasing continuous abatement functions : x f2 — > indexed 
by t G [0,T], we want to find a solution {At}tg[o,T] to 



Af = ttE^ 



T 



1[£t- / rs{As)ds > 



te[o,T]. (9) 



The results of the discrete-time analysis given in [15] suggest that, if the increments 
of the martingale 

{£t:=E'^i£T\Tt)}mo,T] 

are independent and the abatement functions : x f2 — )• are deterministic and 
time independent, then one can reasonably expect that a solution to (Q can have the 
functional form 

At = ait,Xt), te[0,T], 
with an appropriate deterministic function 

a: [0,T] X M ^ M (10) 

and a state process {Xt}fg[o,T] given by 

Xt:=St- [ rs{As)ds, t e [0,T]. (11) 
Jo 

In this section, we demonstrate how this approach enables one to find a solution in 
the framework of diffusion processes. Assume that {VFj}tg[o,T] is a standard Brownian 
motion process (under Q ~ P) and that our {J^t}te[o,T] is the natural filtration of the 
process. In this case, by the martingale representation theorem [16], one must have 

d£t = atdWt 

for some admissible adapted process {crt}t(=io,T]- To ensure that {St}t£io,T] has indepen- 
dent increments, we assume that {at = cr(t)}tg[o,T] is a known deterministic function 
and that we are given continuous non-decreasing time-independent abatement func- 
tions {vf = r}tg[o,T]- To verify the martingale property of the allowance price process 

At = a{t,Xt), te[0,T], 

we use Ito's formula and ( ITT]) to write the stochastic differential of the process as 

dAt = da(t, Xt) 

= 9(i,o)a(t, Xt)dt + 9(o,i)a(t, Xt)dXt + ^9(o,2)«(t, Xt)d[£]t 

= 9(i,o)a(t, Xt)dt - <9(o,i)a(t, Xt)r{a{t, Xt))dt + i9(o,2)«(t, Xt)a'^{t)dt 
+ d(^^i)a{t,Xt)a{t)dWt. 
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Here [£]t stands for the quadratic variation of the martingale {£t}t(^[Q^T] and de- 
notes the respective partial derivatives. Now we observe that the function a is a 
solution to 

9(i,o)tt(t, x) - r{a{t, x))d(^o,i)a{t, x) + ^a'^{t)d^o,2)a{t, x) = (12) 

in (0, T) X M with the boundary condition 

a(T, x) = 7rl(x ^ 0), x G M, (13) 

justified by the digital terminal allowance price, then the thus constructed {Atjtgjo^T] 
will indeed be a martingale that satisfies by construction. Note that a{t, x) satisfies 
the maximum principle, so that ^ a{t, x) ^ vr for all {t, x) G [0, T] x R. For the proof 
of this fact, see Proposition 14. II below. 

The following summarizes the above-presented approach. 

Allowance price in the diffusion framework 

1. Given a continuous non-decreasing function r : [0, oo) — )> [0, oo) and a positive 
function cr(t), t G [0,T] x M, determine a solution a to the boundary value 
problem f|T2|) . f|T3|) . (We assume that a(t) is regular enough to ensure existence 
and uniqueness of the solution.) 

2. Verify that there is a unique strong solution to 

dXt = dSt - r{a{t, Xt))dt, Xo = So. (14) 

3. Introduce the allowance price {At}t£[o,T] by 

At:=a{t,Xt), te[0,T]. 

Having constructed the allowance price process {At}tg[o,T] in this way, one obtains 
a standard procedure for the valuation of European options. Indeed, observe that, 
due to the Markov property of the strong solution to (fT4l) . the fair time t price of a 
European call option written on the allowance price at (maturity) date r G {t, T] is 
given in terms of an appropriate function of the state variable: 

Ct = [{Ar - K)+ I J^t] = [{a{T,Xr) - KT I J-t] = nt,Xt). 

To ensure that {Ct = f^lt, Xt)}t£[o^r] is a martingale, the function : [0, r) x M — > R 
is to be taken as a solution to the linear partial differential equation 



9(i,o)r(t,x) - d^o,i)r{t,x)r{a{t,x)) + ^9(o,2)r(t, x)a2(t) = 0, (15) 
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in (0, r) X M. However, the boundary condition in this case will be 



r{T,x) = {a{T,x)-K) 



(16) 



Summarizing, we obtain the following description for the procedure. 
Valuating a European call in the diffusion framework 

1. Find the function a as above. 

2. Given a strike price K > and maturity time r G [0, T] of a European call, 
calculate as the solution to the boundary problem f[T5]) . ffTBl) . 

3. Given a time t G [0, r] and the allowance price a G [0, tt] at time t, obtain x as 
the solution to a{t,x) = a. 

4. Substitute t and the thus obtained x into the function to obtain the time t 
price of the European call as f"^{t,x). 

Note that one can also estimate directly using Monte Carlo simulations: given a 
strike price K > 0, a. maturity time r G [0, T] of a European call, a time t G [0, r] and 
the allowance price a G [0,7r] at time t, obtain Xt = x a.s the solution to a{t,x) = a. 
Then, using the Markov property of {Xt}, one can evaluate the option price 



by estimating the expectation via simulating a (large enough) number Nmc of copies 
of the random variable X^-. 

Note that closed-form solutions to the non-linear partial differential equation ( !T2|) 
are rarely available. However, a linear abatement function leads to explicit expressions, 
as pointed out in [25] and [31]. We will consider this case as an illustration. 

Example. Given a linear abatement function r(a) = ca, c E (0, oo), and a constant 
diffusion coefficient {at = a}t£[o,T] with a G (0, oo), the partial differential equation 
( !T2|) becomes Burger's equation 



whose solution can be obtained from that of the heat equation via the Hopf-Cole 
transform. 

Namely, a direct calculation shows that if f : [0, T] x M i— R solves the heat equation 



= [{a{T,X,)-K)+ \Xt = x] 



d(i,o)U - cud(o,i)U + ^<9(o,2)M = 0, 



(17) 




(18) 
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then its Hopf-Cole transform 

u := — — (19) 

c V 

solves ( fTTl) . In order to satisfy an original boundary condition 



the boundary value function for ( 1181) must be chosen as 

f (T, x) = 6(a;) := exp | — y j> , a; G R. 

Our digital boundary condition has the form 

u{T, x) = b{x) = 7rl(x ^ 0), x G M, 

so that we should take 

v{T, x) = b{x) = l{x < 0) + l{x ^ 0)6""^'=/'"', xeR. 

Denoting by N(/i, s^) the normal distribution with mean /i and variance and by 
and f{y) = ^'{y) the standard normal distribution function and its density, 
respectively, we can write: 

v{t,x) = j b{y)N{x,a\T-t)){dy) 

/O poo 
N{x,a'{T-t))idy)+ / e-^-^/-'N(x, ^^(T - 
oo Jo 

^ f -X \ fx-cTx{T-t)\ r CTT 7r2c2(T-t)l , , 
= ^(-l^]+^\ h^P o,. ^ • (20) 



The derivative of this expression with respect to x is equal to 



(9(o,i)f (t, x) = , H , exp < -X + 



Next we use f ll9p to calculate 



a{t, x) = ^ ' / ^ \ \ (t, x) G [0, T] X M. (22) 
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Figure 1: The functions a(t, ■) displayed for t = 1.9, 1.6, 1.3, 1.0, 0.7, 0.4 and calculated 
using ([2U])-(I22D with T = 2, a = 4, c = 0.02 and vr = 100. 



The shape of this function is depicted in Figure [H 

To illustrate our valuation procedure, consider the following parameters: time to 
compliance date T = 2, diffusion coefficient cr = 4, penalty vr = 100 and a linear 
abatement function r(a) = ca with c = 0.02. Taking t = 0, we consider a family of 
European calls with the same strike price K = 25, but different maturity times r. 
Suppose that the initial allowance price is equal to the strike price a = Aq = 25 (the 
so-called at-the-money situation), which is attained by a{0,x) = 25 with x ~ —2.434. 
Next, we determine the call prices C(0,r) at time t = for different maturity times 
r G [0,T]. Independently of the model, the price of an expiring call with r = must 
be equal to zero, whereas the longest-maturity call with r = T must have the price 
AqItt — K)/7r = 25 X 0.75 = 18.75. (Note that, because of the digital payoff, such a 
call is equivalent to 0.75 allowances). That is, the call prices increase with contract's 
maturity from to 18.75. The shape of this curve is obtained using a crude Monte 
Carlo procedure. We have simulated Nmc = 10^ i.i.d. copies of the random variable 
Xr- We used the forward Euler method to generate i.i.d. copies of {Xt}, with time 
step 0.02. The curve is presented in Figure [2], together with 95% confidence intervals. 



4 Allowance price in the jump-diffusion setting 

In this section, we describe a situation where the market can suddenly change due to 
the presence of jumps. We propose a framework suitable to the present practice of the 
EU ETS where the member states negotiate their allowance allocations and the market 
needs to adapt to the new situation. For instance, revised decision on the amount of 
allocated certificates yields a jump of allowance market prices. Also, a sudden change 
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Figure 2: The prices of European calls depending on their maturity times. 



in demand and/or price for fuel can result in the pollution levels changing dramatically 
which impacts on the allowance prices. 

We discuss how to solve problem ([8]) in the jump-diffusion setting using a similar 
approach to the one employed in Section [31 Namely, we choose a candidate for iQ of 
the form 

At = a{t,Xt), 

where a and Xf are as in fllUp and flTTl) . respectively. But now we assume that the 
martingale 

is modeled using a general jump-diffusion process adapted to a filtration {J^t}t£[o,T]- 
In this setting, the process {Xt}t&[o,T] is given by the following stochastic differential 
equation: 

dXt = -r{a{t,Xt))dt + a{t,Xt)dWt + / a{t-, Xt.,y)ip, - q,)idy x dt), (23) 

where Mq = M \ {0}, {W^t}te[o,r] is a Brownian motion adapted to the filtration 
{Tt}te[o,T], and p^, is an independent of {W^t}te[o,r] {-^t}te[o,r]-adapted random Pois- 
son measure with intensity qu{dy x dt) = \iy{dy)dt, v being a probability distribution 
on M and A G (0, oo) a positive constant. Expression (12^ allows us to model jumps 
with great flexibility, as the distribution of the jump is state and time dependent. It is 
known that, under suitable Lipschitz and growth conditions on r, a and a, stochastic 
differential equation fl23|) possesses a unique strong solution [H]. 

We again use Ito's lemma for At = a{t,Xt) to find conditions on a under which 
the process {At}te[o,T] will be a martingale under Q. Let t = T — t. For convenience. 
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we introduce a new function /3(r, x) = a{T — t,x) and, instead of ffT^ . obtain the 
following nonlinear partial integro-differential equation for /3: 

d{ifi)P{T,x) = -r(/3(r,x))(9(o,i)/3(r,a;) + ^cr^(r, a;)9(o,2)/3(r, x) 

+-^y" [f^i-r^x + a{T,x,y)) - f3{T,x) - a{T,x,y)d(^o^i)f3{T,x)]u{dy), (24) 

(r, x) G [0, T] X M, with the boundary condition 

/3(0,x) = 7rl(x ^ 0), xeM. (25) 

Under certain assumptions, jump-diffusion models for option pricing lead to partial 
integro-differential equations, see e.g. [UIH]- Equation (I24p differs from them in that it 
has a nonlinear coefficient r(/3(r, x)). We assume that there exists a classical solution 
(3 e Ci'2([0,T) X R) to the problem (EH), ([25]). It is not difficuh to see that it satisfies 
the maximum principle. 

Proposition 4.1. Solution /3(r, x) to equation [24\ ) with initial boundary condition 
/3(0,x) = h{x) satisfies the maximum principle: 

inf h{z) ^ /3(r, x) ^ sup h{z), (r, x) G [0, T] x M. (26) 

2 2 

Proof. Since with a(t, x) = /3(T — t,x), where /3 satisfies 1^^, the process = a{t, Xt) 
is a Q-martingale. Also, the process {Xt} given by f l23|) is a Markov process. Thus, 
one has a.s. 

a{t,Xt) = E{At I Tt) = E{a{T,XT) \ Ft) = E(a(T,Xr) | Xt) 

and so, due to regularity of a, for r = T — t, 

/3(r,x) = a(t,x) = E(a(T,XT) \ Xt = x) 

= E(/3(0, Xt) I = x) = E(/i(Xt) | X* = x), 

which immediately implies ( 126|) . □ 

Note that Proposition 14.11 also proves the maximum principle for solution to the 
boundary value problem f lT2|) . f|T3|) . 

We will see that, under a mild condition, the maximum principle remains valid for 
the discretized problem as well (see Proposition 14.21 below) . 

Analytical solution to partial integro-differential equations can only be obtained in 
a few special cases, so in most situations one can only solve them numerically. In this 
section, we discuss a discretization of (124|) using the finite difference method and prove 
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that the discretized equation has a unique solution. Finally, we illustrate the approach 
by providing a numerical example. 

First we need to truncate the domain of x. Let := {x e M : |x| < /}, where the 
bound / > of the domain can be chosen so that the probability that the process {Xt} 
leaves Di during the time interval [0, T] given it starts at Xq does not exceed a given 
small e > 0. This procedure will be illustrated at the end of the section. 

Next, we restrict the domain of integration for the integral term on the RHS of 
flM|) to an interval [Ki , , chosen such that the error made due to the truncation 
also remains small. For a good choice of the terminals Kj, one can refer to the study 
presented in [3]. 

The numerical solution will be computed on a discrete grid. Let denote the total 
number of discrete a:-values and M the total number of r-values we want to use in 
the grid for the numerical solution, so that the step sizes in x and r are respectively 

= 21/N and = T/M. Put Xi = -l + iA^, Tn = nA^, for i e Z and n = 0, . . . , M. 
We use /3f = /3{Tn,Xi) for the values of /3 on this grid. Because of the presence 
of the non local term on the RHS of equation (l24l) . one needs to define /3 outside 
[0, T] X Di. We chose the simplest and most intuitive approach and set /3f := g{xi) for 
i ^ {0, . . . ,iV - 1}, where 



9[x) 



71 if X ^ I, 
if X ^ -/. 



The partial derivatives are replaced by the respective finite differences: 



9(1,0) 13 {Tn,Xi) 
9(0,1) (3 {Tn,Xi) 
9(l,0)(3{Tn,Xi) 



A, ' 



Following an approach similar to the one used in [TO], we consider the same step size 
Aj. to approximate the integral term, and choose Ji and J2 such that [Ki,K2] C 
[(Ji — l/2)Aa;, (J2 + l/2)Aj.], which leads to the approximation 

^ J2 

A / /3(r„, Xi + a(r„, Xi, y))v{dy) ~ A ^ ^jl3]*(i,j,n) 



where 



0+l/2)A, 

Uj:= j v{dy) 

(j-l/2)A, 
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and 

Similarly, one has 



j*{ij,n):= arg min + a(r„, Xj, x^)) - (A;A^ - 

k 



J2 

where a^j := a{Tn,Xi,Xj). We set := ^ arj'^j convenience. 

So now we are looking for a solution {/?"} on the grid such that, for n = 0, . . . , M—1, 

= 2^ (-n+i,a:.) ^-^^ r(A ) ^- 

+ ^ E - A/??-^^ - A^^^— ^sr^ (27) 



with 



/3° = 7rl(xi ^ 0) for i e Z, /3f = ^(xi) for i < and i > N. 

Let S* := minS" and a*^ := min (T^(r„, Xj), where the minima are taken over all 
i e {0, . . . , iV - 1} and n G {0, . . . , M}. 

Proposition 4.2. T/ie discretized problem has a unique solution {/3"}. //, in addition, 
cr^ is bounded away from zero and the discrete grid is such that 

-S*A, ^ — , 
2A ' 

then the solution satisfies the maximum principle: 

< ^ TT for anyieZ andne{0,...,M- 1}, 

where tt is the penalty per unit of pollutant not covered by the initial allocation. 

This proposition ensures that the discretized problem is well posed. In particular, 
it shows that the allowance price process is always positive and does not exceed the 
level of penalty fixed by the regulator at any time during the compliance period. 

Proof. We follow the steps in the proof from [10]. For n G {0, . . . , M — 1} and i G 
{0,...,iV — 1}, equation ( j271) can be rewritten as 

-F{p-)A^/3i\' + (1 + Gnp?)K)pr' - ^rAr/3r+Y (28) 

= /3r + AA.^z/,/3]l(,,^.„), 
j=Jl 
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where 

G-(6-) - + + + A (30) 

= 2(A.)^ • ^''^ 
This is a hnear system and can be written as 

M(n)/3"+^ = y„, for all n = 0, . . . , M - 1, (32) 
where = (/3o"+\ • • • , P%t\) G M^, Yn = (y^', • • • , G with components 

Jr 

= /3r + AA. ^^-/^rC..,"), ^ = 0,...,N-2 
j=Ji 

and 

Jr 

given by the right hand side of fl28l) . Given /3", the matrix M(n) G M^^^ is tridiagonal: 
the elements on its main diagonal are the terms (l + G"(/3f )A^), the elements of the fist 
diagonal above it are given by —if "A,- and the elements of the fist diagonal below the 
main diagonal are given by — F"(/3")At-. Furthermore, M(?t,) is diagonally dominant, 
which can be seen from the relation 



G'^{f3^) = FJ'{(3^) + H^ + X for z = 0, . . . , - 1 and n = 0, . . . M 

and the non- negativity of the coefficients (l29l) - (!3T]) . Therefore, given /J", the linear sys- 
tem (!32|) possesses a unique solution /3"'^^, see [18]. Hence the existence and uniqueness 
of the solution to (!27|) follows by induction. 

Now we will use induction in n to show that the maximum principle holds. We will 
only prove that (3f are non-negative, as the argument can be easily adapted to prove 
that the values remain bounded by tt. 

We want to show that, for any A,- > and A^. > 0, 

/3f ^0 for all i e Z and n E {0, . . . , M - 1}. 

For n = 0, this is obvious from the shape of the boundary condition f|25p and the 
definition of g. For induction step, assume that /3f ^ for all i E 1^, but there exists 
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an Zo G Zi such that (3'^^'^ < 0. By definition of g, io E {0, . . . , N — 1} since we would 



have Z?"'*'^ = g{xif^) ^ otherwise. One can choose iq such that 



= min < 0. 

Under assumption (14.21) . F", G" and iff are all non-negative. It follows from G^{I3^) = 
F^{I3^) + H^ + \ that 

^ -FpiPDA^Pzti + (1 + Gmo)Ar)f3:;' - HrA^p:;^\ 

j=Ji 

which is a contradiction. The proposition is proved. □ 

In conclusion of this section, we give a method for choosing the domain boundary / 
and provide a numerical illustration. Let us focus on a special case where a(r, x,y) = y 
holds for all ?/ G M, which corresponds to a compensated compound Poisson process. 
Under this assumption, {Xi}tg[o,T] follows a jump-diffusion process given by 

ft pt Nt 

X, = Xo- r{ais,X,))ds+ / X,)rfiy, + V - AE(ri)t, 

where {Nt}te[o,oo[ stands for a Poisson process with intensity A and {yjjjgN is a sequence 
of independent identically distributed random variables following the same u = N{0, 1), 
all the components {VTt}, {Nf} and {Yj} of the model being independent of each other. 
Denote the compensated jump component by 

Nt Nt 

and the martingale part by Mt := a{s, Xs)dWs + Jt- Suppose for simplicity that 
Xo = 0. 

The reduction function r is non-decreasing and the non-negative function ■) is 
bounded by vr since it satisfies the maximum principle. Thus, for the drift term we 
have ^ 

0^ / r{a{s,Xs))ds ^r{-K)t. (33) 



Next, we will use Kolmogorov-Doob inequality to bound the martingale part of {^t}, 
for which we need a bound for the second moment of Mt. Observe that the covariance 
between the diffusion and jump terms is given by 



^[jt a{s,Xs)dw}j =E Jt^i^j^ a{s,X,)dW,\{Tj,Y,}^ 



(34) 
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where Tj is the time of the j-th jump Yj. Considering the inner expectation given 
{Tj = tj,Yj = Hj}, one can easily see that it is equal to zero, and hence the expected 
value fl34l) is zero as well. It follows that 



Clearly ^{Jf) = XEYit. Next, assume that a satisfies the standard linear growth 
condition 



Then 



a'^{s, x) ^ a + bx'^ , a, 6 ^ 0, xG 



E (^J^ a{s, Xs)dWs^ = E(t2(s, Xs)ds ^ E {a + hXl) ds 



t 

2, 



= at + b I EXtds, 

and so we have 







EM^ ^ XEY^H + at + b f EX'^ds. 

Jo 



Since we see from fl33|) that 

Mt - r(7r)t ^Xt^Mt (35) 

and EMf = 0, one has 

EX^ ^ r\7r)f + EM^ ^ {XEY^ + a)t + r'^^Ti)^ + b [ EX^ds. 

Jo 

Hence Gronwall's lemma (see (2.10) and (2.11), Chapter 5 in [TB]) applied to {Xt} 
provides the upper bound 

EX^ ^ ^ (« + ^^Yi + e'* - ^ (« + ^^n' + 2r'(7r)t + =: Kf (36) 

The Kolmogorov-Doob inequality applied to the submartingale {—Xt} together with 
(|36|) yields 

P ( inf Xt^-l] = P ( sup (-Xt) ^ / K =: ei 

for all / > 0. Noting that 

EMf ^ Kt + r^iTT)t^ =: Ct, 
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Figure 3: Function a(t,-) for t = 0.2,0.4,0.6,0.8,1 calculated by discretization ( 128|) 
with parameters given in the text. 




Figure 4: A typical path of the jump- diffusion process {Xj}tg[o,T]- 

one obtains similarly from the right inequality in f l35|) that 

P ( sup ^ / ] ^ P ( sup Mt^l) ^ Ct^' =: 62, 

which provides a way to determine the domain of truncation [— /, /] for our problem for 
given accuracy levels ei, €2- 

Suppose that the time unit is one year and assume that the time to compliance 
is T = 1. Set the penalty level at vr = 1 and take for simplicity a constant function 
a{t,x) = a in the diffusion term. We consider cr = 1 here. The compound Poisson 
process is realized with intensity rate A = 1 and a standard normal jump distribution 
ly = iV(0, 1). Finally, we suppose that the cumulative reduction function r is linear: 
r(a) = a for all a G [0,7r]. In this case, for ei = €2 = 0.05, we obtain the domain 
truncation boundary / ^ 11. 

We consider a sequence of discretization schemes with the time step = 0.02 
which gives M = 50 time points on the grid. Choosing the same space discretization 
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Allowance price 




Figure 5: The corresponding realization of the allowance price process {At = 
a{t, Xt)}te[Q,T] for the {Xt}te[o,T] obtained in Figure H]. 



Aj. = 0.02, we obtain = 2//0.02 = 100/ space points, depending on the truncation 
[—1, 1] of the space region. To control the error made by the truncation, the Li-norm 
of the difference between two solutions corresponding to the truncations 

at levels ±/ and ±/, respectively, is considered: 

M 

EEiA"(0-/3r(oi, i<l 

i n=0 



where the summation is taken over all grid points on a fixed sub-region [0, T] x [— rf, d] 
with d < I < I. Numerical experiments show that, for d = 11 and / = 20, / = 30, 
the error is of order 10"^. We therefore conclude that the solution calculated 
for the space region [—30, 30] is accurate. Figure [3] shows the shape of the function 
a obtained in this way. The graph in Figure H] depicts a typical realization of the 
process {Xt}fg[o,T], calculated using the forward Euler method [T7]. Figure E] shows the 
corresponding realization of the price process {At = Xt)}(g[o,T]- 

Figure E] displays the function a at time t = 0.8 obtained for various volatility levels 
a and jump rates A. It turns out that, when the volatility is small, the price process 
tends faster to its boundary value: agents do not expect big changes in pollution emis- 
sions and therefore in the allowance price. If there is a shortage (excess) in allowance 
credits before maturity, the market expects that there will also be a shortage (excess) 
at t = T. In the presence of jumps or a large diffusion coefficient, the allowance price 
tends to converge slower as one can expect a sudden increase or decrease in pollution 
emissions which impacts directly on the allowance price. 
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Figure 6: Function a(0.8, ■) for different values of a and A. 

5 Conclusions 

Tfie growing evidence of the cost of climate change justifies the introduction of large- 
scale measures. Emission trading schemes become increasingly important. The pro- 
posed operation of newly designed cap-and-trade mechanisms within a multi-period 
setting and with inter-connection to other markets adds a notable complexity to such 
systems, and so quantitative understanding of emission trading schemes becomes in- 
creasingly challenging. For instance, the problems of market design, the emission reduc- 
tion performance, the optimization of allowance allocation procedures, the individual 
risk management and the valuation of emission-related financial instruments need to be 
addressed within a sound mathematical framework, which we aim to approach in this 
work. Based on results from equilibrium analysis, we focus on the simplest situation 
of a one-period market to show how the risk-neutral evolution of emission allowances 
can be described in terms of jump-diffusion processes. Although the resulting partial 
integro differential equations are non-linear, we provide numerically stable and fast val- 
uation procedures which yield reliable numerical schemes for valuation of derivatives 
of the fast-growing family of emission-related financial assets. 
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